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Abstract: Bacterial and fungal communities play critical roles in reestablishing vegetation structure, 
function and biodiversity in ecosystem restoration in arid and semi-arid areas. However, the long-term 
successional changes in bacterial and fungal communities that occur with artificial vegetation development 
are not fully understood. In this study, we investigated the successional changes in bacterial and fungal 
communities in Caragana korshinskii Kom. plantation over a period of 50 a (6, 12, 18, 40 and 50 a) and 
their relationships with key soil environmental factors in a typical agro-pastoral ecotone, northern China. 
The results showed that bacterial and fungal diversities (x- and B-diversity) were significantly affected by 
plantation age; moreover, the change in fungal community was more evident than that in bacterial 
community. Soil samples from 12 a plantation had the highest (P<0.05) bacterial and fungal «-diversity (i.e., 
abundance-based coverage estimator (ACE) and Chaol index) at 0-10 cm depth compared with other 
samples. However, soil samples from plantation at the late recovery stage (40-50 a) had the highest 
«-diversity at 10-20 cm depth. Soil bacterial community was not significantly affected by plantation age at 
the genus level; but, soil fungal community was significantly affected at the genus level. Overall, Mortierella 
and Chaetomium were the dominant genera at natural recovery stage (0 a); Inocybe was the dominant genus at 
the early recovery stage (6—12 a); Inocybe and Mortierella were the dominant genera at the mid-recovery 
stage (12—40 a); And Mortierella, Cladosporium and Humicola were the dominant genera at the late recovery 
stage (40-50 a). Redundancy analysis (RDA) showed that 6-glucosidase activity, total nitrogen and soil 
organic carbon were closely associated with bacterial community composition, while alkaline phosphatase, 
urease activity and total nitrogen were associated with fungal community composition, indicating that 
changes in enzyme activity and soil nutrients were the most important determinants of dominant genera. 
Furthermore, pathogenic microorganisms (Cladosporium and Humicola) were dominant in soils from 40—50 
a plantation, which may affect plant growth, resulting in the decline of C. korshinskii plantation. Overall, 
the findings of this study improve the understanding of ecological patterns of bacterial and fungal 
communities in artificial vegetation and provide an important scientific basis for comprehensive ecological 
restoration management in arid and semi-arid areas. 
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1 Introduction 


Microbes play critical roles in soil nutrient cycle, structural formation and plant interaction. These 
roles are important for reestablishing soil microbial function and biodiversity during ecosystem 
restoration (Jim, 2009). For instance, soil microbes can promote the formation of humic acid that 
is an important component for soil improvement and aboveground plant growth (Schinel, 1995; 
Harris, 2003; Sun et al., 2017). Moreover, some microbes can help improve stress resistance and 
resource efficiency in terrestrial plants (Mendes et al., 2011), while soil-borne pathogenic 
microbes can affect the health of plants (Wall et al., 2015). Generally, microbes are more sensitive 
to external stress than animals and plants (Panikov, 1999); therefore, external factors have a 
considerable impact on soil microbial community structure and function. Soil microbial 
community composition and diversity could be a reflection of changes in soil environment and 
plant community during the process of vegetation restoration (Yeates, 1979; Wall and Moore, 1999; 
Wu et al., 2008; Hu and Guo, 2012; Frouz et al., 2016; Du et al., 2018). In arid and semi-arid 
areas, soil microbes play an important role in ecosystem protection and in maintaining stable 
productivity of plant community. 

Since the past decades, the area lost to desertification has been increasing globally (He et al., 
2005; Zhao et al., 2018). China is one of the most severely affected countries by land 
desertification (Ma and Zhou, 2007). Chinese government has implemented a series of ecological 
management policies for ecological restoration, such as returning farmland for forestry, wind and 
sand source control of Beijing and Tianjin, and Three-North Shelter Forest Construction (Ma et 
al., 2006). 

Caragana korshinskii Kom. belongs to the Leguminosae family. This plant with an extremely 
developed root system is highly resistant to drought, heat, cold, salinity and alkaline conditions 
and presents a xeromorphic structure. Furthermore, C. korshinskii plays an important role in water 
and soil conservation, vegetation restoration and ecological environment improvement. C. 
korshinskii is the main shrub species naturally growing in harsh environments (Zhang et al., 2014), 
and it has been planted on a large scale in the desert areas of Northwest China (Niu, 2003). As a 
result of efforts to reduce desertification, the area of artificial vegetation in China is the largest 
and accounts for approximately one-third of the world's afforestation area (Zhang and Gao, 2000). 
The large-scale construction of artificial vegetation could influence soil microbial structure and 
soil property. For instance, Wang et al. (2013) found that Pinus sylvestris var. mongolica Litv. and 
Populus simonii Carr plantations significantly improved soil physical-chemical property and 
increased microbiological activity; consequently, higher soil microbial abundance, microbial 
biomass carbon and enzymatic activity were recorded in these plantations compared with mobile 
dunes in Horgin Sandy Land. Similarly, Liu et al. (2009) reported that long-term planting of 
Cunninghamia lanceolata could affect soil microbial community by altering soil property, which 
could consequently influence nutrient cycle in forest plantation. Additionally, long-term 
cultivation of Pinus massoniana (25 a) can improve soil microbial community stability and 
microbial function (Zhao et al., 2020). However, the impact of artificial vegetation on the 
ecological environment is broad and an important subject to study. Presently, although some 
researchers have investigated microorganisms in C. korshinskii plantation, few studies have 
directly explored bacterial and fungal changes over a long period. 

The agro-pastoral ecotone in northwestern Shanxi Province is one of the most vulnerable 
ecosystems in northern China because of intensive human activity and a fragile natural 
environment, which is susceptible to severe soil erosion and land degradation. Additionally, this 
area has complex and diverse land use types, comprising of agricultural lands (i.e., corn, potato 
and bean), artificial vegetation (i.e., C. korshinskii, Pinus tabuliformis, Populus and Salix 
matsudana Koidz) and grassland. Therefore, owing to the various environmental gradients (i.e., 
time series, spatial heterogeneity and land use types) in agro-pastoral ecotone, this is ideal site for 
studying the succession of microbial community and their relationships with key environmental 
factors. In this study, we investigated the changes of soil bacterial and fungal communities of C. 
korshinskii plantation over a 50-a period (6, 12, 18, 40 and 50 a) and their relationships with key 
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soil environmental factors in a typical agro-pastoral ecotone in northern China. The result might 
provide a scientific basis for the management and ecological construction of C. korshinskii 
plantation in arid and semi-arid areas. 


2 Material and methods 


2.1 Study area 


The study was conducted in Shizuitou Village, Wuzhai County, Xinzhou City, Shanxi Province, 
China (38°44'—-39°17’N, 111°28’-113°00'E). The average altitude was approximately 1400 m a.s.l. 
This area has a temperate continental climate, with annual precipitation of 450-500 mm, annual 
evaporation of 1784 mm and annual average wind speed of 3 m/s. January is the coldest (-13.3°C) 
and July is the hottest (20.1°C). The number of windy days in spring is more than 36 d. The soil 
type is classified as loessial soil, which is characterized by a loose texture and low fertility. 
Natural vegetation in this area was originally a grass-dominated steppe, but during the past 50 a, 
natural vegetation had been substantially altered due to prolonged heavy grazing by livestock. 
Reestablishment of artificial vegetation on degraded land is one of the most effective ecological 
practices in this region. Since the 1970s, artificial C. korshinskii plantation has been established in 
this area. Moreover, the plantation was established in different time periods, making the area ideal 
for studying changes in soil and vegetation. We divided artificial vegetation into four types 
according to the plantation age: natural recovery (0 a, CK), early stage (6-12 a), mid-stage (12—40 
a) and late stage (40-50 a). Changes in morphological characteristics of C. korshinskii plantation 
at different ages are shown in Table 1. 


Table 1 Changes in morphological characteristics and root biomass of Caragana korshinskii plantation at 
different ages 


Plantation age 


Index 
6a 12a 18a 40a 50a 
Crown (m°) 1.29+0.14° 2.12+0.16° 4.55+0.25° 8.22+0.52° 6.56+0.22° 
Plant height (cm) 111.56+15.66° 131.35416.85% 178.47+37.20° 216.00+39.64" 202.50+46.97” 
Root bi : 
PARE 0.28+0.02° 0.32+0.01°4 0.34+0.03% 0.36+0.01° 0.36+0.01° 
(kg/plant) 


Note: different lowercase letters within the same row indicate significant differences among different age groups (P<0.05 level, Duncan 
tests). Means+SD. 


2.2 Sampling site and soil sampling 


Soil samples were taken from two depths (0—10 and 10-20 cm) under the shrub canopy at six 
different ages (i.e., 0, 6, 12, 18, 40 and 50 a) from July to October in 2019. A hand auger boring 
was used for sampling. The sample plots were located on the flat land between hills, and they had 
similar site conditions (i.e., soil matrix, vegetation and climate). Natural recovery land was 
selected as control (0 a). We used a nested experimental design to randomly select three sampling 
sites with a similar slope degree (4°—8°) and aspect (South) of each plantation age. The slope 
positions of sampling sites were located top-slope, mid-slope and down-slope of each plantation 
age, which were separated by at least 200 m. Then, three 20 mx20 m quadrats, separated by at 
least 25 m, were set at each sampling site. At each quadrat, a 3-m-wide and 20-m-long transect 
was set up, and five soil samples under C. korshinskii and between C. korshinskii were collected 
along the transect and then pooled to form two composite samples for each quadrat. For each age, 
18 soil samples were collected, and 108 soil samples were totally obtained. Each sample was 
divided into two sub-samples: one was stored at —20°C for soil DNA extraction and enzyme 
activity analysis, and the other was air-dried for soil physical-chemical analysis. Soil 
B-glucosidase, alkaline phosphatase (ALP) and urease activity were measured as described by Liu 
et al. (2019). 


2.3 Soil physical-chemical analyses 


Air-dried soil samples were used for soil physical-chemical analyses. Soil moisture, soil salt 
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content and pH were measured using a soil five-parameter analyzer (COMBI 5000, Hamburg, 
Germany). Soil nutrient content including soil organic carbon (SOC) and total nitrogen (TN) were 
measured by the method of Xi et al. (2015). 


2.4 Soil deoxyribonucleic acid (DNA) extraction, polymerase chain reaction (PCR) and 
Miseq sequencing 

We extracted total genomic DNA from 0.5 g samples using Power Soil DNA Isolation Kit (MO 
BIO Laboratories, Carlsbad, CA, USA) according to the manufacturer's protocol. To profile soil 
bacterial community, we amplified the V3—V4 region of the 16S rRNA (ribosomal ribonucleic 
acid) gene with the forward primer 5'-ACTCCTACGGGAGGCAGCA-3' and reverse primer 
5'-GGACTACHVGGGTWTCTAAT-3’, combined with adapter sequences and barcode sequences. 
For fungal community, we amplified the ITS—ITS1 region of the fungal 18S rRNA region using 
the forward primer 5'-CTTGGTCATTTAGAGGAAGTAA-3' and reverse primer 
5'-GCTGCGTTCTTCATCGATGC-3'. PCR amplification was performed in a total volume of 50 
uL, which contained 10 uL buffer, 0.2 uL Q5 high-fidelity DNA polymerase, 10 uL high GC 
enhancer, 1 uL dNTP (deoxy-ribonucleoside triphosphate), 10 uM of each primer and 60 ng 
genomic DNA. Thermal cycling conditions were as follows: an initial denaturation at 95°C for 5 
min, followed by 15 cycles at 95°C for 1 min, 50°C for 1 min and 72°C for 1 min, with a final 
extension at 72°C for 7 min and storage at 4°C. The amplified products were subjected to 2% 
agarose gel electrophoresis, and the amplified products were subjected to [lumina MiSeq 
high-throughput sequencing and analysis. The sequencing and bioinformatics services used in this 
study were completed by the Biomarker Technologies Corporation, Beijing, China. 


2.5 Data analysis 


Poor quality sequences, including sequences <50 bp in length and bases with quality <20, were 
screened out using QUME (Quantitative Insights Into Microbial Ecology) software version 1.9.1, 
and chimeric sequences were identified and removed using UCHIME 8.1 software to obtain a 
high-quality tag sequences. The qualified sequences were clustered into OTUs (operational 
taxonomic units) at a similarity of 97% using USEARCH 10.0. Furthermore, species annotation 
and taxonomic analysis of OTUs were performed using Silva (http://www.arb-silva.de/) and Unite 
(Release 8.0, https://unite.ut.ee/) databases for bacteria and fungi, respectively. 

The a-diversity measurements of soil bacteria and fungi, including the number of observed 
OTUs, Chaol, ACE (abundance-based coverage estimator) and Shannon and Simpson diversity 
indices were determined using QIIME software. We used Adonis analyses based on the binary 
Jaccard distance matrixes of bacterial and fungal communities to determine the significance of 
microbial community differences among different soil groups using the vegan package in R 
software version 4.1. 

One-way analysis of variance (ANOVA) and Spearman's correlation were performed using 
SPSS software v13.0, and the significance level was set to P<0.05 level. Canoco v4.5 was used to 
perform RDA (Redundancy analysis) on the relationship of bacterial and fungal communities with 
soil properties. Graphical illustrations were prepared using OriginPro 2017 software. 


3 Results 


3.1 Changes in bacterial and fungal compositions 

3.1.1 Bacterial composition 

At the phylum level, plantation age did not significantly affect the dominant bacterial phyla at 
0-10 and 10-20 cm depths (P>0.05). The dominant phyla were Proteobacteria, Acidobacteria, 
Actinobacteria, Gemmatimonadetes and Chloroflexi, with Proteobacteria accounting for more 
than 25.12%—40.09% of total reads in samples. Bacteroidetes, Rokubacteria, Planctomycetes, 
Nitrospirae and other bacteria constituted a minor proportion (Fig. 1). 
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Fig. 1 Relative abundance of dominant bacteria (a and b) and fungi (c and d) at the phylum level. CK, control. 


At the class level, plantation age did not significantly affect the dominant classes in bacterial 
community (P>0.05). Alphaproteobacteria, Gammaproteobacteria, Subgroup _6, 
Blastocatellia Subgroup_4, Gemmatimonadetes and Acidimicrobiia accounted for major 
proportions in samples at different plantation ages (Fig. 2). 
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Fig.2 Relative abundance of dominant bacteria (a and b) and fungi (c and d) at the class level. CK, control. 
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At the OTUs level, plantation age did not significantly affect the composition of dominant 
OTUs (P>0.05). The dominant OTUs were affiliated with Sphingomonadaceae, uncultured_ 
bacterium c Subgroup 6 and Gemmatimonadaceae, followed by some OTUs within 
Pyrinomonadaceae, Nitrosomonadaceae and uncultured_batetiun_o IMCC26256 (Table 2). At 
0-10 cm depth, Sphingomonadaceae had the highest relative abundance, followed by 
uncultured_bacterium_c_ Subgroup 6 and Gemmatimonadaceae. At 10-20 cm depth, 
uncultured_bactertum_c Subgroup 6 had the highest relative abundance followed by 
Gemmatimonadaceae and Sphingomonadaceae. 


Table 2 Dominant bacteria operational taxonomic units (OTUs) and their relative abundance at different 
plantation ages 


Relative abundance (% 
OTUs Affiliation Depth a 
(cm) CK 6a 12a 18a 40a 50a 
f Sphingomonadales, 0-10 18.141.4  19.543.3 13.342.0 14.641.4 12.8414 16.1413 
Sphingomonadaceae 10-20 14.043.6 7.92.1 4,140.7 5.440.6 6.2+1.0 7.10.9 
Uncultured_bacteritum 0-10 9.0+0.4 7.9+0.4 10.0+0.3 9.3+0.5 9.740.3 10.7+0.5 
2 _c Subgroup 6, 

‘Acidobacteria 10-20 10.0+0.2 8.90.2 10.3+0.4  9.440.7  — 10.0+0.3  10.9+0.2 
, Gemmatimonadales, 0-10 7.8+0.3 6.6+0.8 7.9+0.5 7.5+0.1 7840.3  6.8+0.2 
Gemmatimonadaceae 10-20 8.8+0.4 6.7+1.0 9.5+0.0 8.9+40.1 9.440.4 8.8+0.4 
j Blastocatellales, 0-10 5.1+0.0 5.3+0.2 5.4+0.6 4.8+0.5 5440.3 5,440.2 
Pyrinomonadaceae 10-20 4540.4 4.0+0.1 4.140.3 3.2+0.1 3.940.2  4.9+0.3 
: Nitrosomonadales: 0-10 2.4+0.1 2.6+0.3 3.140.3 3.140.4 3.140.3 2.4+0.2 
Nitrosomonadaceae 10-20 2.9+0.5 3.6+£0.2 4.2+0.1 4340.2 4340.1 3.30.1 
Actinobacteria, 0-10 2.6+0.2 2.50.1 2.80.2 2.90.3 2.70.1 2.10.1 

6 Uncultured_batetiun 
Actinobacteria, 0-10 2.3+0.4 2.3+0.1 2.740.5 2.3+0.2 2.60.3 1.50.1 

7 Uncultured_batetiun 
c MB-A2-108 10-20 3.5+0.6 4.7+0.3 5.6+0.3 5.7+0.4 5.2+0.4 4,240.4 


Note: CK, control; Mean+SE. 

At the genus level, plantation age did not significantly affect the composition of dominant 
genera (P>0.05). The dominant bacterial genera at 0-10 cm depth were 
uncultured_bacterium_c_ Subgroup _6, uncultured_bacterium_f_Gemmatimonadaceae, 
Sphingomonas and RB41. The dominant bacteria genera at 10-20 cm depth were 
uncultured_bacterium_c_Subgroup_6, uncultured_bacterium_f_Gemmatimonadaceae, 
Sphingomonas and uncultured_bactertum_c_MB-A2-108 (Fig. 3). 


3.1.2 Fungal composition 


At the phylum level, plantation age significantly affected the composition of dominant phyla in 
fungal community (P<0.05). At natural recovery (CK), the early (6-12 a) and late (40-50 a) 
stages, Ascomycota and Mortierellomycota were dominant phyla at 0-10 and 10-20 cm depths, 
and their relative abundances were 47.13%—64.11% and 7.62%-—18.05%, respectively. However, 
at the mid-stage (12—40 a), Ascomycota and Basidiomycota were dominant phyla, and their 
relative abundances were 46.96%—72.56% and 16.08%-—24.75%, respectively. Glomeromycota, 
Chytridiomycota, Zoopagomycota, Olpidiomycota, Kickxellomycota and unclassified fungi 
constituted a minor proportion (Fig. 1). 

At the class level, plantation age significantly affected the composition and relative abundances 
of dominant classes in fungal community (P<0.05). At natural recovery stage (0 a), 
Dothideomycetes and Sordariomycetes were dominant classes at 0-10 and 10-20 cm depths, 
respectively. At the early and mid-stages (6-12 and 12—40 a), Agaricomycetes and Pezizomycetes 
were dominant classes. At the late stage (40-50 a), Sordariomycetes, Leotiomycetes and 
Mortierellomycetes were the dominant classes (Fig. 2). 

At the OTUs level, plantation age significantly affected the composition of dominant OTUs 


ChinaX iv ERAT 


MA Gailing et al.: Succession of soil bacterial and fungal communities of ... 


202111.00045v1 


naXiv 


(P<0.05). At natural recovery stage (CK), Nectriaceae, Chaetomiaceae and Mortierellaceae were 
dominant at 0-10 and 10-20 cm depths. At the early stage (6-12 a), Chaetomiaceae, Inocybaceae 
and Pleosporales_fam Incertaesedis were dominant. At the mid-stage (12—40 a), Mortierellaceae, 
Inocybaceae and Thelephoraceae were dominant. At the late stage (40-50 a), Mortierellaceae, 
Chaetomiaceae, Cladosporiaceae and Nectriaceae were dominant (Table 3). 
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Fig.3 Relative abundance of dominant bacteria (a and b) and fungi (c and d) at the genus level. CK, control. 


Table 3 Dominant fungi operational taxonomic units (OTUs) and their relative abundance at different plantation 


ages 
Relative abund: % 
OUTs Affiliation Depth oes 
(cm) CK 6a 12a 18a 40a 50a 
0-10 11.240.6 17.2+7.1 6.61.4 10.3+5.6 8.2+2.5 10.843.1 
1 Sordariales, Chaetomiaceae 
10-20 15.5+2.6 3.6+0.8 FIELS 6.72.8 10.943.1 17.6+5.9 
0-10 7.2£1.3 2.6£0.7 12.6+3.0 13.1+5.9 VI 742.7 17.141.3 
2 Mortierellales, Mortierellaceae 
10-20 = 16.9+8.9 1.8+0.1 10.9+2.6 5.4+2.0 11.8+3.9 9.71.6 
0-10 2.31.2 9.42.9 6.8+5.3 11.6+6.5 11.146.3 0.76.3 
3 Agaricales, Inocybaceae 
10-20 5.2+2.5  37.0+17.9 10.9+410.0 12.9410.2 24.0+13.0 0.80.3 
0-10 5.5+2.1 0.7+0.2 5.741.4 7.041.3 6.62.0 8.641.8 
4 Hypocreales, Nectriaceae 
10-20 3.0+0.1 1.340.3 4.7£1.7 3.10.9 4.50.9 6.11.3 
P Pleosporales, 0-10 0.20.0 23.7+13.6 0.70.4 0.30.2 0.10.0 2.5£1.5 
Pleosporales_fam Incertaesedis 10-20 0.140.0 30.9417.8 0.10.1 0.1+0.0 0.340.2 0.20.1 
7 Cladosporiales, 0-10 1.90.0 4.443.7 6.1+4.4 3.0+40.6 3.8+40.9 5.2413 
Cladosporiaceae 10-20 1.640.2 0.80.3 2.2+0.8 2.140.3 3.0£0.8  7.8+0.9 
0-10 0.10.0 2.3+0.6 4,142.7 5.44+4.2 3.6+1.6 0.4+0.3 
7 Thelephorales, Thelephoraceae 
10-20 0.3+0.1 2.4+0.1 F252 5.43.1 5.62.8 0.40.3 


Note: CK, control; Mean+SE. 
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At the genus level, plantation age significantly affected the dominant fungal genera in the soils 
(P<0.05). At natural recovery stage (CK), Mortierella and Chaetomium were dominant genera at 
0-10 and 10-20 cm depths, whereas Jnocybe was dominant genus at the early recovery stage 
(6-12 a). At the mid-recovery stage (12—40 a), Inocybe and Mortierella were dominant genera; 
Mortierella, Cladosporium and Humicola were dominant genera at the late recovery stage (40-50 
a) (Fig. 3). 

3.2 Bacterial and fungal a-diversity 


Plantation age significantly affected (P<0.05) a-diversity of bacterial community at 0-10 and 
10—20 cm depths. For 0-10 cm depth, the ACE and Chaol indices of soil from 12 a plantation 
were significantly higher than those of other samples. However, the Simpson and Shannon indices 
were not significantly affected by plantation age. For 10-20 cm depth, the ACE and Chaol 
indices of soil from 50 a plantation were significantly higher than those of other samples. 
Additionally, the Simpson and Shannon indices of soils collected from 6 a plantation were higher 
than those of other samples (Fig. 4). 
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Fig. 4 Variations in ACE (al—a4), Chaol (b1—b4), Simpson (cl—c4) and Shannon (d1—d4) indices of bacterial 
and fungal communities. * and ** indicate significant difference among plantation ages at P<0.05 and P<0.01 
levels, respectively. CK, control. 


Plantation age significantly affected a-diversity of soil fungal community (P<0.05). For 0-10 
cm depth, the ACE, Chaol, Simpson and Shannon indices of soils from 12 a plantation were 
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significantly higher than those of other samples. For 10-20 cm depth, the ACE and Shannon 
indices of soils from 50 a plantation were significantly higher than those of other samples. 
However, the Chaol index of soil from 40 a plantation was significantly higher than that of other 
samples. The highest Simpson index was obtained in soils from 12 a plantation (Fig. 4). Overall, 
fungal community was more affected by plantation age than bacterial community. However, 
bacterial community was richer and more diverse. 


3.3 Bacterial and fungal B-diversity 


Beta diversity analysis is usually used to compare the similarity of different groups in species 
diversity. The binary Jaccard algorithm calculates the distance between groups, which is 
expressed as B-diversity. The higher the B-diversity, the greater the species difference between 
groups. The B-diversity of bacterial and fungal communities was significantly affected (P<0.05) 
by plantation age, and it exhibited a quadratic pattern. At 0-10 and 10-20 cm depths, the 
difference in bacterial community between groups is greater than that of within groups (7=0.202 
and 0.273, respectively). The differences in bacterial community between groups peaked at 12 a 
and 40 a plantations, respectively. At 0-10 and 10-20 cm depths, the difference in fungal 
community between groups is greater than that of within groups (7=0.505 and 0.585, 
respectively). The differences in fungal community between groups peaked at 18 a and 12 a 
plantations, respectively (Fig. 5). 
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Fig. 5 B-diversity (binary Jaccard) of bacterial (a and b) and fungal (c and d) communities. "All between" 
represents the beta distance data of all samples between groups. CK, control. 


3.4 Changes in soil properties 


3.4.1 Soil physical-chemical properties 


Soil physical-chemical properties, such as soil moisture, total salt contents and pH were 
significantly affected by plantation age. Soil moisture exhibited a quadratic trend, i.e., it peaked in 
soils from 12 a plantation and decreased afterwards at both 0—10 and 10-20 cm depths. Similarly, 
soil salt content showed a quadratic trend, i.e., it peaked in soils from 12 a plantation and 
decreased afterwards. At 0-10 cm soil depth, soil pH showed a decreasing trend, and the lowest 
and highest pH values were obtained in samples from 12 and 50 a plantations, respectively. At 
10-20 cm soil depth, soil pH was not significantly affected by plantation age (Fig. 6). 
Furthermore, SOC, TN and C:N ratio were significantly affected by plantation age. At 0-10 cm 
depth, the lowest and highest SOC values were obtained in soils from 6 and 40 a plantations, 
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respectively, whereas SOC at 10-20 cm depth was not significantly affected by plantation age. At 
0-10 cm depth, the lowest and highest TN values were obtained in soils from 6 and 50 a 
plantations, respectively. However, at 10-20 cm depth, the lowest and highest TN values were 
obtained in soils from 40 and 50 a plantations. At both 0-10 and 10-20 cm depths, the lowest and 
highest C:N ratios were obtained in soils from 6 a and 50 a plantations, respectively (Fig. 6). 
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Fig. 6 Changes in soil physical-chemical properties (a, soil moisture; b, soil salt content; c, pH value; d, soil 
organic carbon; e, total nitrogen; f, C:N ratio) of Caragana korshinskii plantation. Bars are standard errors. 
Different lowercase letters indicate significant differences among different plantation ages at P<0.05 level. CK, 
control. 


3.4.2 Soil biological properties 

Soil -glucosidase and ALP activities were significantly affected by plantation age. At 0-10 cm 
depth, the lowest and highest B-glucosidase activities were obtained in CK soil and that from 12 a 
plantation, respectively. At 10—20 cm depth, the lowest and highest B-glucosidase activities were 
obtained in soils from the 6 a and 12 a plantations, respectively. ALP exhibited a quadratic trend. 
The lowest value was obtained in soils from 18 a plantation, and increased afterwards. Soil urease 
activity was not significantly affected by plantation age at 0-10 and 10-20 cm soil depths (Fig. 7). 
Overall, the activity of soil B-glucosidase was higher than those of ALP and urease. 


3.5 Relationships of compositions of bacterial and fungal communities with soil properties 


Physical (soil moisture), chemical (SOC, TN and soil salt content) and biological properties of 
soil (ALP, B-glucosidase and urease) at different plantation ages were used as main environmental 
factors in RDA analysis. 

For bacterial community, 57.9% of cumulative variation could be explained by RDA analysis. 
The first and second axes explained 64.1% and 91.5% of cumulative variations, respectively. SOC, 
TN, B-glucosidase and soil salt content were closely related to the dominant genus composition of 
bacterial community. 

For fungal community, 45.6% of cumulative variation was explained by RDA analysis. The 
first and second axes explained 88.3% and 96.4% of cumulative variations, respectively. The 
genera Mortierella and Chaetomium were dominant at natural recovery stage (CK). Soil moisture 
and B-glucosidase activity were positively correlated with Chaetomium abundance. Urease 
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Fig. 7 Changes in soil enzyme activities (a, B-glucosidase; b, alkaline phosphatase; c, urease) of Caragana 
korshinskii plantation. Bars are standard errors. Different lowercase letters indicate significant differences among 
different plantation ages at P<0.05 level. CK, control. 


activity and TN were positively correlated with Mortierella spp. TN and urease activity were 
positively correlated with dominant species of fungal community (Mortierella, Cladosporium and 
Humicola) at the late recovery stage (40-50 a), whereas soil salt content, SOC and ALP were 
positively correlated with dominant genus of fungal community (Inocybe and Mortierella) at the 
early and mid-recovery stages (6-12 and 12—40 a) (Fig. 8). 
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Fig. 8 Redundancy analysis (RDA) plots showing the relationships of bacterial (a) and fungal (b) communities 
with soil properties. SOC, soil organic carbon; TN, total nitrogen; ALP; alkaline phosphatase. 


4 Discussion 


This study revealed the successional patterns of bacterial and fungal communities over a period of 
50 a in C. korshinskii plantation, located in northwestern Shanxi Province, China. Compared with 
bacterial community, fungal community showed more significant changes in dominant species 
and diversity over time. Changes in fungal community could be attributed to changes in soil 
properties. Particularly, soil nutrients (i.e., TN and SOC), enzyme activity (ALP) and plantation 
age were highly associated with fungal community succession. 


4.1 Changes in dominant species of fungal and bacterial communities 


This study showed that dominant fungal taxa were significantly affected by plantation age 
(P<0.05). At natural recovery stage, the dominant classes were mainly microorganisms with a 
high tolerance to stress (i.e., Sordariomycetes and Dothideomycetes), and at the early and 
mid-recovery stages, the dominant species were mainly microorganisms with a moderate 
tolerance to stress, which can encompass harsh to moderately rich soil environments (i.e., 
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Agaricomycetes and Pezizomycetes). However, at the later recovery stage, the dominant species 
had a low tolerance to stress (i.e., Leotiomycetes and Mortierellomycetes). The dominant species 
found in our study area are consistent with those reported in other arid and semi-arid regions 
(Bastida et al., 2014; Maestre et al., 2015; Martirosyan et al., 2016; Rao et al., 2016). These 
results demonstrate that in arid and semi-arid areas, the soil environment is relatively single, so 
that the diversity of soil microorganisms is not high and the species composition is relatively 
simple. At the late recovery stages, soil microorganisms had a low tolerance to stress, so that the 
abundance of beneficial microorganisms decreased, while that of pathological microorganisms 
increased, resulting in the degradation of C. korshinskii plantation. 

Indicator of fungal genera was also identified at different plantation ages. The dominant fungal 
genera in our study area were Mortierella, Chaetomium (natural recovery stage), Inocybe (early 
and mid-recovery stage), Cladosporium and Humicola (late recovery stage). At the late recovery 
stage, beneficial microorganisms (Inocybe) were replaced by pathological microorganisms 
(Cladosporium and Humicola). These significant changes have likely taken place in response to 
alterations in the soil's physical and chemical properties. At the late recovery stage, soil moisture 
content decreased and salt content increased. Drought and salt stress decreases enzyme activity 
and soil nutrients, with corresponding changes of the composition of dominant microbial species. 
At the genus level, dominant fungal genera were different from those present in arid and 
semi-arid areas. These differences may reflect habitat heterogeneity and variations in topography 
(i.e., top-slope and down-slope). Under different land uses, microbial community respond 
differently to environmental factors, and the number and composition of soil microbes changes 
accordingly (Drenovsky et al., 2010). Moreover, the changes in dominant species at different 
plantation ages indicate that fungal community succession occurs largely through the adjustments 
of dominant species, and only species better adapted to the soil environment at a specific stage 
could exist. 

For bacterial community, the dominant species were not significantly affected by plantation age 
(P>0.05), which is consistent with those from other arid and semi-arid areas (Neilson et al., 2012; 
Taketani et al., 2015; Yao et al., 2017). The bacterial compositions at the genus level were very 
different from that reported in other areas. For example, Arthrobacter, Bacteroides, 
Faecalibacterium, Sphingomonas and Gaiella were reported as dominant in northern Ningxia, 
China (Wang et al., 2021); while Sphingomonas, Microbacterium, Bradyrhizobium and 
Pedomicrobium were dominant bacteria in the middle Loess Plateau, China (Wang et al., 2019). 
These differences may reflect the habitat heterogeneity, including climate, soil and vegetation. 
The loess hilly area has complex habitats and a strong spatial heterogeneity. However, the 
composition of bacterial community is not significantly affected by the increase of plantation age, 
reflecting that the structure of soil bacterial community is relatively stable. 


4.2 Changes in bacterial and fungal diversities 


At 0-10 cm depth, a-diversity of bacterial and fungal communities from 12 a plantation were 
significantly higher than those from other ages, but at 10-20 cm depth, a-diversity from 40-50 a 
plantation was opposite, which suggested that richness and diversity indices of bacteria and fungi 
appeared lagging behind as the soil depth increased. Similar result was confirmed by Liu et al. 
(2019). The bacterial and fungal communities in deep soil might use environmental resources in 
more restrictive ways than those in shallow soil. 

Indices of a- and B-diversity of bacterial and fungal communities showed similar patterns at 
different plantation ages. Indices of a- and B-diversity of bacterial and fungal communities tended 
to increase at early stage and then decreased, at later stages, which was similar to those of other 
studies in arid and semi-arid regions (Shi et al., 2014; Tedersoo et al., 2014; Liu et al., 2021). For 
example, Hunt et al. (2003) found that species richness of soil microorganisms in North America 
gradually increased during the early stages of plantation but decreased as the forest aged. Higher 
biomass and litter formation provide better conditions for the accumulation of soil nutrients at 
early plantation stage. With an increase in soil nutrients and an improved soil structure, soil 
microorganisms become more abundant. However, as succession progresses, soil moisture 
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decreases and precipitation events can no longer replenish soil moisture in deep soil, resulting in a 
decrease in the number and diversity of microorganisms (Yan and Cao, 2010). Therefore, we 
suggested that artificial C. korshinskii should be stumped after 20 a plantation. Studies have 
shown that a reasonable stumping time for C. korshinskii is 11-15 a (Cheng et al., 2009). 
However, an appropriate stumping time should consider the local topography, climate and other 
environmental factors (Li et al., 2014). 


4.3 Key environmental factors associated with bacterial and fungal communities 


This study indicated that enzyme activities were correlated with dominant bacterial and fungal 
species at different plantation ages. For bacterial community, B-glucosidase, SOC and TN were 
positively correlated with bacterial community. For fungal community, B-glucosidase, soil 
moisture, ALP and soil salt content were positively correlated with fungal community. Studies 
showed that most bacteria prefer habitats with rich soil nutrients but less salinity, while most 
fungi tend to exhibit an opposite trend (Wei et al., 2020). Our study confirmed that SOC and TN 
might play important roles as microbial substrates. The effects of environmental factors on soil 
microbiome are dynamic and these factors do not act in isolation. For instance, soil bacterial 
B-diversity first increased and then decreased with increasing N input (Liu et al., 2021). 
Additionally, the changes in dominant species, microbial diversity and species relative 
abundances could be partly attributed to osmotic and ion toxic effects, resulting from increasing 
soil salt content that may affect soil moisture availability and enzyme activity (Yan et al., 2015). 
In present study, the change of fungal community diversity and species abundance were higher 
than those of bacterial community, indicating that bacterial community was stable and fungi 
might play important roles with the increase in plantation ages. 

In addition to soil environmental factors, biotic interaction in bacterial and fungal communities 
is also important (Pontarp and Petchey, 2016; Ning et al., 2019). At early and middle stages, 
Inocybe can participate in decomposition of organic matter, promote circulation of C, N and P, 
and can also degrade a variety of environmentally harmful substances. At late stage, Mortierella, 
Cladosporium and Humicola replaced Inocybe as dominant genus. These fungi may affect plant 
growth, resulting in the decline of C. korshinskii because Cladosporium and Humicola can cause 
plant diseases such as plant leaf spots, leaf mold and stem rot. Several studies have shown that 
continuous cropping are mainly associated with the increase in crop diseases and subsequent 
decrease in the yield (Hu et al., 2006; Liu et al., 2009). Therefore, future management of C. 
korshinskii in late recovery stage should focus on the suppression of pathological 
microorganisms. 

This study showed that the fungi:bacteria ratio increased with increase in plantation age (Fig. 
9). Similar result was confirmed by Helm et al. (1996) and Zhao et al. (2011). Most bacteria 
prefer habitats rich in soil nutrients but with less salinity, whereas most fungi tend to exhibit the 
opposite trend. Fungi have a higher tolerance to water stress than bacteria (Wang et al., 2019). 
Under harsh environmental conditions, fungi have stronger vitality than bacteria (Chen et al., 
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Fig.9 Changes of fungi:bacteria ratio at different plantation ages 
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2007). Therefore, the soil fertility declined with the continuous plantation of C. korshinskii and 
fungi became the dominant microorganisms. 


5 Conclusions 


Our study revealed the succession of soil bacterial and fungal communities in C. korshinskii 
plantation over a period of 50 a. This vegetation succession was reflected in compositional 
changes, a- and B-diversity of bacterial and fungal communities. Enzyme activity and soil 
nutrients were highly correlated with the composition of bacterial and fungal communities. At the 
late recovery stage (40-50 a), soils were dominated by pathogenic microorganisms 
(Cladosporium and Humicola) with potentially adverse effects on plant growth, resulting in the 
degradation of C. korshinskii plantation. Thus, future strategies for forest recovery and 
establishment should focus on improving soil's physical and chemical environment through 
artificial intervention or introducing beneficial microbial species that have antagonistic effects on 
those pathogenic microorganisms. 
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